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Abstract 

We derive the nonperturbative effect in type OB string theory, which is defined by 
taking the double scaling limit of a one-matrix model with a two-cut eigenvalue 
distribution. However, the string equation thus derived cannot determine the 
nonperturbative effect completely, at least without specifying unknown boundary 
CN ■ conditions. The nonperturbative contribution to the free energy comes from in- 

stantons in such models. We determine by direct computation in the matrix model 
an overall factor of the instanton contribution, which cannot be determined by the 
string equation itself. We prove that it is universal in the sense that it is inde- 
pendent of the detailed structure of potentials in the matrix model. It turns out 
to be a purely imaginary number and therefore can be interpreted as a quantity 
related to instability of the D-brane in type string theory. We also comment on 
a relation between our result and boundary conditions for the string equation. 
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1 Introduction 



The nonperturbative effect in string theory can be studied using matrix models. In 
particular, the noncritical string theory, which is a simplified model of string theory, is 
exactly solvable via matrix models l a . The string equation, which can be derived from 
them, contains the nonperturbative effect of the noncritical string theory |2j. On the 
other hand, study of the Liouville theory [3] enables us to find the effect of the D-brane, 
which is the nonperturbative effect of string theory and can be identified with the effect 
that appears in the string equation. In [4. , we have shown that the string equation does 
not describe the nonperturbative effect completely, at least in c = noncritical string 
theory. To obtain the whole nonperturbative effect, it is necessary to study the matrix 
model directly. 

Recently, a matrix model that corresponds to the noncritical string with worldsheet 
supersymmetry has been proposed Ej. For the c = noncritical string, which is 
described as two-dimensional pure supergravity on the worldsheet [7], we consider the 
double scaling limit around the Gross- Witten transition |Sj. This critical point can 
be found in the one-matrix model with a two-cut eigenvalue distribution. In string 
theoretical interpretation, the one-matrix model with two cuts corresponds to the NSR 
string theory of type OB. This matrix model can be solved with the string equation. 
However, as in the case of c = string theory, it does not contain the nonperturbative 
effect completely. 

In this paper we study the nonperturbative effect of type OB string theory by analyz- 
ing the matrix model with a two-cut eigenvalue distribution. We compute the effect of 
instantons directly in the matrix model, which corresponds to the D-brane in the string 
theory. The result is summarized as follows. 

From the string equation, the nonperturbative effect in the free energy is obtained 

as 
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where C cannot be determined from the string equation by itself without specifying 
unknown boundary conditions. From the direct computation using the matrix model, 
we can determine the constant C as 

Moreover, it is shown that this value is universal, namely, it does not depend on the 
detailed structure in the potential of the matrix model. Because it is purely imaginary, 
this nonperturbative effect is related to the instability of the D-brane. 
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The paper is organized as follows. In section we identify the instantons in the 
matrix model and the contribution from instantons to the free energy. In section El we 
compute the effect of the instantons using the method of orthogonal polynomials. In 
section we take the double scaling limit and consider the universal behavior of the 
effect of instantons. In section we present the conclusions. The appendices show the 
details of calculations. 

2 Instanton in one-matrix model 

In this section, we consider a one-matrix model. We discuss how an instanton contributes 
to the partition function and the free energy. The one-matrix model with a one-cut 
eigenvalue distribution corresponds to c = noncritical string theory, while that with 
a two-cut distribution corresponds to c = type OB string theory [Hj. In both cases, 
an instanton can be interpreted as the ZZ-brane In the one-cut case, the instanton 
can be described as a configuration in which all eigenvalues are at the minimum of the 
effective potential except that a single eigenvalue is at its maximum. This description 
can be extended to the case of two-cut distribution. In this case, the effective potential 
has two separated minima. Half of the eigenvalues are in one of these minima and the 
other half are in the other minimum except that a single eigenvalue is at the maximum. 
In the one-matrix model, the partition function is given by 



Here, $ is an iV x iV Hermitian matrix. Hereafter we consider the case where the 
potential V(x) is invariant under x — > —x and thus the eigenvalue distribution of <3> has 
this Z2 symmetry. 

Diagonalizing the matrix <3>, the partition function can be expressed as 



Here, A(A) = Yli<j — ^j) * s ^ ne Vandermonde determinant. We concentrate on the 
iV-th eigenvalue A at, and represent it as x. The other N — 1 eigenvalues can be regarded 
as those of a (N — 1) x (N — 1) matrix. The partition function of an N x N matrix 




(2.1) 




(2.2) 
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model can be expressed using an (N — 1) x (N — 1) matrix model as 

N-l 



/r " ' —N ) V(Xi)-NV(x) 

dx / n dA, A^A) ]J(x- A,) 2 e fei (2.3a) 

= y da; (det(x - $) 2 > 7V _ 1 e~ NV ^ (2.3b) 

= Z N _ ± J dxe~ NV ^ x \ (2.3c) 

Here, the subscript "JV — 1" indicates that the quantities concerned are those in the 
(N — 1) x (TV — 1) matrix model, and the expectation value (O) is defined as 

(O) = ^ J d$Oe- NtrV W. (2.4) 

In the large-iV limit, the system of an (N — 1) x (TV — 1) matrix is the same as the 
system of an N x N matrix. Hence, we can use the standard N x N matrix model to 
calculate these expectation values. 

The effective potential V e g(x) defined above can be expressed in terms of connected 
diagrams. After some algebra, we obtain 

NV cS (x) = NV(x) - 2 (trlog(:r - $)> c - ± ((2trlog(x - $)) 2 ) c - • • • . (2.5) 

Here, the subscript "c" indicates the connected part. In the large- N limit, the first and 
second terms are of order N and third term is of order N°. If we restrict ourselves to 
the leading order of N, the terms other than the first two can be neglected. Using the 
resolvent 1 

R(x) = (^tr— ^ = 1 (y'(x) - y/V»(x)+p(x)) , (2.6) 
the equation can be expressed as follows: 

NV cS (x) = NV^ ] + Vg + ^ + • • • , (2.7a) 
= V(x) - 1 (2trlog(x - $)) (2.7b) 

/X px 
dx'R(x') = Re / dx'^V'^x') +p(x). (2.7c) 

The resolvent R(x) has the cut on the real axis. If x is on the cut, the effective potential 
becomes constant and the eigenvalue density p(x) takes a nonzero value. Interpreting 
this in physical terms, we are considering the iV-th eigenvalue x and the other N — 1 



1 Here the branch of the square root is chosen so that R(x) ~l/iasi^oo. 
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eigenvalues are those of $ whose distribution is expressed by p(x). In the cut, where the 
N — 1 eigenvalues are distributed, the forces from them acting on the iV-th eigenvalue 
cancel each other. The effective potential V e f[(x) has is at the minimum over the entire 
cut. From the standpoint of the original system of the N x N matrix including the 
iV-th eigenvalue, it is natural that the eigenvalue density p(x) should not change due 
to the iV-th eigenvalue at the leading order of N; that is, the back reaction is at the 
subleading order. Hence, in the integration with respect to x in (j23J), most of the 
contribution comes from the region where x is inside the cut. In the case of one cut, the 
integration over this region gives the partition function of the N x N matrix system. 
To extend this to type OB string theory, we should consider the case of two cuts. In 
the case of two cuts, there are two minima of the effective potential. Because the 
potential under consideration is Z 2 symmetric, these two minima should be symmetric 
under x — > —x and give the same contribution to the partition function. Hence, we 
can deal with these two regions together as "inside the cut." There is another nonzero 
contribution from the region where x lies outside the cut. Because there is a maximum 
of the effective potential, we should take this into account. If the iV-th eigenvalue x 
lies outside the cut, the eigenvalue density of the N x N matrix system differs from 
that of the (N — 1) x (JV — 1) matrix system. It can be identified as the "instanton" 
characterized by the configuration where x is located at the local maximum of the 
effective potential. The resolvent is related to the disk amplitude of the Liouville field 
theory with the boundary condition corresponding to the FZZT-brane ^0]- The effect 
of the instanton comes from the (integrated) resolvent with some fixed value of the 
cosmological constant corresponding to the local maximum of the effective potential. It 
is related to the ZZ-brane jH]- 

Now, we divide the integration region into two parts, namely, inside the cut and 
outside the cut, as 

Z = Z N „! I dx (det(x - $) 2 ) e~ NV{x) 

J inside the cut 

+ Z N _ 1 f dx(det(x-$) 2 )e- NV( > x \ (2.8) 

J outside the cut 

The second term comes from the integration outside the cut, and is identified as the 
contribution of the instanton. So far, we have restricted ourselves to the iV-th eigenvalue 
x and identified the instanton as the configuration where x lies outside the cut. However, 
there are N eigenvalues and all other N — 1 eigenvalues can be possibly outside the cut 
as well. Hence, there is an n-instanton sector; that is, n eigenvalues lie outside the cut. 
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The partition function can be expressed the sum of those in the n-instanton sector as 

Z = ^(0-inst) + ^(1-inst) + ^(2-inst) + . . . ? ( 2 .9a) 

-N^2v(X t ) 



^(o-inst) = /" JJdA i A 2 (A)e" 

J inside the cut 



= zflT ] / dx <det(x - $) 2 ) ((Mnst) e~ NV(x \ (2.9b) 

«/ inside the cut 

Z (l-inst) = ^(0-inst) f ^ ^ et(x _ $)2 ^(0-inst) ^ 

J outside the cut 

Here, the superscript "(n-inst)" indicates the n-instanton sector. In the partition func- 
tion ^( - mst ) and the expectation value ^(9}^ 0_mst ^ in the O-instanton sector, all eigenvalues 
are inside the cut. Hereafter we will omit the superscript "(0-inst)" in the expectation 
value. The factor N in front of the 1-instanton sector partition function ()2.9c|) reflects 
the number of ways of specifying an eigenvalue that lies outside the cut. If we con- 
sider all n-instanton sectors, neglecting interaction between instantons, which is valid 
in the large- N limit, the partition function can be expressed in terms of those of the 
O-instanton and the 1-instanton sector 



e F 



(y(l-inst) \ 

Z {1 -"" 1> f,A„. ,di(det(a;-4>) 2 )e~ m ' W 

]\T J outside the cut \ V ' ' Icy i ni \ 

P - Z (o-inst) T dx (det(x - $) 2 > e-^W ' 1 ' 

J inside the cut N v / ' 

The additional term // is the contribution from instantons. We regard this term as the 
"chemical potential" of the instanton. This effect corresponds to the ZZ-brane in the 
noncritical string theory. 

3 Effective potential and orthogonal polynomials 

In the previous section we have seen that the partition function has contributions which 
come from the multi-instanton sectors, that are characterized as configurations where 
some of the eigenvalues are located at the maximum of the effective potential. On 
the other hand, in order to compute the chemical potential of the instanton, which is 
nothing but the instanton effect in the free energy, it is sufficient to consider the effect of 
only one instanton. This amounts to considering (det(x — $) 2 ) in the O-instanton sector 
and performing the integration with respect to x as shown in (J2.10bj) . However, the 
integration of (det(x — $) 2 ) inside the cut generally contains divergent contributions in 
the subleading order of N 4|. Our expectation here is that these divergences cancel the 
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overall N in (|2.10b|) to make the chemical potential finite. In order to confirm this, we 
have to retain the iV-dependence in the computation. For this purpose, it is appropriate 
to use the method of orthogonal polynomials, because it is available for any N. 

We begin with definitions and properties of the orthogonal polynomials. The parti- 
tion function of the one-matrix model can be expressed in terms of orthogonal polyno- 
mial P n (x) as 

Z = / TT dA * detP n (A n /)detP m (A m Oe < . (3.1) 

/ nn' mm' 

J i 

Here, the orthogonal polynomial P n {x) = x n + 0{x n ~ 1 ) satisfies the orthogonality con- 
dition 

(P n (x),P m (x)) = J dxP n (x)P m (x)e- NV{x) = hj nm . (3.2) 

Using this inner product of the orthogonal polynomials, the partition function can be 
expressed as 

N-l 

Z = N\ det(P n , P m ) =N\h£T[ f»~ n . (3.3) 

nm -*- A 

n=l 

Here, f n = h n /h n -\. The orthogonal polynomials can be determined by recursion 
relations 

xP n {x) = X nm P m (x) = P n+1 {x) + s n P n {x) + r„P n _i(x), (3.4a) 

P' n { X ) = VnmPm(x) = [NV (X)} nm P rn (x) . (3.4b) 

It can be easily seen that r n = f n . Eliminating r n and s n , we will obtain a differential 
equation that determines P n (x). However, we use these relations in a slightly different 
way. Because the free energy is expressed in terms of r n 's, we should determine them. 
This can be done using (J3.4bj) . Indeed, r n and s n can be determined by (|3.4bJ) . then 
using ()3.4a|) . P n (x) can be expressed in terms of r n and s n . 

In the large-^ limit, it is natural that any quantity f n with index n can be approx- 
imated by a continuous function /(£) with £ = n/N. In fact, /„ becomes continuous 
in the case of one cut. In order to do this, however, the values of f n and f n+ \ should 
become closer in the large- iV limit. In the case of two cuts, /„ cannot be approximated 
by a continuous function, but if we consider f n with even n or odd n separately, they 
can be approximated by different continuous functions [TT]. We use f n to indicate that 
the index n is even, and f n for odd n. In the large- N limit, they are approximated by 
different functions as ^ 

f = l^ = L^ ( n: 6Ven ' ) (3 5) 

\fn = /(0 (n : odd) ' 1 • ] 
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In this limit, a summation over the index n can be approximated by an integration. If 
we want to calculate up to the next-to- leading order in 1/N, we should use the Euler- 
Mclaurin summation formula. For example, for even N 

A N f 1+ ^ ~ N f 1 ~ 

£/» = ¥ / , d ^) + T / d ^^' ( 3 - 6 ) 

n=0 J°~7f ^° 

In this manner we treat corrections of the next-to-leading order in a systematic way. 
Now, we compute the chemical potential of the instanton 

f ... t dx(det(x-$) 2 )e- w W . 

r .j . t ,dx(det(x - $)2) e -^W ' r 

J inside the cut > \ / / 

Here, it is easy to show that the expectation value (det(x — $)) can be identified with 

Pn 

P n (x) = (det(ar-$)) n , (3.7) 

where the subscript "n" again indicates a quantity for an n x n matrix system. Because 
the quantity under consideration is not the expectation value of the trace, but of the 
determinant, the large- N factorization does not hold in this case. In fact, if we define 
D n (x) as 

D n (x) = (det(x-$) 2 ) n , (3.8) 

it satisfies a recursion relation 

D n = P*(x)+r n D n _ 1 . (3.9) 

Using this relation recursively, we obtain 

D N (x) = P 2 N {x) + tnPI^x) + ■ ■ ■ + r N ■ ■ ■ rxP 2 (x). (3.10) 

This formula enables us to evaluate the chemical potential. Substituting this formula 
into (|2.10b[) . the chemical potential can be expressed in terms of the orthogonal poly- 
nomials, which can be determined by (|3.4j) . and we will obtain the definite value of the 
chemical potential. 

When we evaluate (det(x — $) 2 ), the result will be different depending on whether 
x is inside the cut or outside the cut. This difference can be described as follows. The 
orthogonal polynomial P n can be expressed as the expectation value in the nxn matrix 
system as in (|3.7)1 . and thus if x is inside the cut of this system, P n (x) is oscillating 
rapidly, otherwise P n (x) is monotonic as x n . As n increases, the oscillatory region 
becomes wider. 
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When x is inside the cut, we define n*(x) as the minimum value of n such that 
x is inside this oscillatory region of P n (x). The normalized orthogonal polynomials 
i>n{x) — Pn(x)/Vhn take values of the same order for n > n*(x), while those for 
n < n*(x) ip n {x) become small so that the contributions from the last n*(x) terms 
containing P n (x) (n < n*(x)) in (j3.10|) are damped exponentially. Hence, we can neglect 
the term of P n with n < n*(x) in ()3.1()|) . On the other hand, if x is outside the cut, P n (x) 
for all n are not oscillating but monotonic with respect to x as x n and we cannot neglect 
the latter terms with P n (x) for n < n*(x). Thus, in the computation of (|3.10j) . we 
should consider these two cases, namely, inside the cut and outside the cut, separately. 

First, we consider the case in which x is outside the cut. In this case, the largest 
contribution to (det(x — $) 2 ) comes from the first term in (j3.10|) and contributions from 
latter terms become smaller like a geometric series. Using the ratio of the orthogonal 
polynomials 

k n _ Pn{%) 



Pn-l(x)' 



D N = (det(x — $) 2 ) can be expressed as 



D N = P 2 N {x) [I + r N e~ 2kN + r^-ie" 2 ^- 2 ^- 1 + ■ 



]■ 



(3.11) 



(3.12) 



In the case of one cut, we can approximate this by a geometric series 

1 



D 



N 



E^w^^ 2 *1"(i+o(4)) 

n=0 ^ \ / / 



(3.13) 



In the case of two cuts, however, we should distinguish a quantity with even n and odd 
n. Hence, we approximate this up to 0(1/N) by a sum of two geometric series as 



for even N, and 



D 



N 



?1=0 



N ( 1 + r N e 



r N r N e 



-2(k N +k N ) 



P 



N 



X 



1 + r N e' 



-2k* 



1 — rjyrN exp 



-2(k N + k 



(3.14a) 
(3.14b) 



D 



N 



71=0 



2 N ( 1 + r N e~^ 



r N r N e 



-2(k N +k N )] V 



P 



A' 



X 



1 + r N e' 



-2k N 



r N r N exp 



-2(k N + k 



N 



(3.15a) 
(3.15b) 
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for odd N. In the case of two cuts, we obtain (see Appendix |XJ) 



1.(0) 4- —k^ + 
N 



7.(0) 



[x 2 + r r , 



r n ) 2 - 4f n r n | 



2x 

e fc ™ ' = — | (x 2 - f n + f n ) ± a/ (x 2 - f n - r n ) 2 - 4r n r r 

^ + k ( n ] = log [(X 2 - f n - F n ) 2 - 4f n f n ] , 

where the double sign is understood to be + for x 2 > r n + r n + 2\Jr n r n so that e 
as | a; | — > oo, and — for a; 2 < r n + r „ — 2y/r n r n so that e fc ' <l) does not diverge at x = 
Using these relations, we obtain 

(x 2 =f fjy ± Fat) exp 



(3.16a) 
(3.16b) 
(3.16c) 
(3.16d) 

',(0) 



(0) 



(x 2 — fjv — tat) 2 — irjsrTN 



1 + 0| iv 



d£ U (0) (0 + ^ (0) (0 ) = 2 / dx'ii!(x). 



(3.17a) 



(3.17b) 



Here, the sign "±" distinguishes the cases of even N and odd N. The relation between 
W(x) and resolvent i?(x) can be checked with a concrete potential V(x), or in the double 
scaling limit. We will examine this point later. 

Next, we consider the case where x is inside the cut. In this case, the contribution 
from P n (x) for n < n*(x) is exponentially suppressed. Hence, = (det(x — $) 2 ) can 
be approximated as 



D 



N 



N-n-1 

(•r) J | r N . 

m=0 



(3.18) 



n=n* (x) 

In the large- N limit, the summation over n can be replaced by an integration. By 
definition, x is in the oscillatory region of P n {x) with n > n*{x). Because its frequency 
is of the order N as shown in Appendix El in order to calculate P n (x) 2 in the large- N 
limit it is sufficient to take the average of this oscillation, which amounts to dividing 
the amplitude of P n (x) 2 by 2. In the case of two cuts, generally, we should distinguish a 
quantity for even n and a quantity for odd n. However, after we average this oscillation, 
i/j n (x) = P n (x) I \fh^ can be approximated by a continuous function of £ = n/N whether 
n is even or odd. Finally we should take the extra factor 2 into account. This factor 
comes from a relative normalization between the P n {x) inside the cut and P n (x) outside 
the cut as in the continuity formula in the usual WKB approximation [4. Using the 
asymptotic formula for if) n (x) inside the cut given in Appendix El (|3.18|) takes a value 
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up to 0(1/N) in the exponent 

N N-n-l 



D N = 22 h N^l(x) I [ r N _ m (3.19a) 

n=n, m=0 

= 2e NReW ^^N I df X = (3.19b) 



V 4 ^(0-(^ 2 -^0-^(0) 5 

27rNp(x)^e N{v{x)+Wo \ (3.19c) 



where = n*/N and is again related to the resolvent as in (J3.17bj) . Because the 

real part of the resolvent becomes V'(x) inside the cut, ReW(x) = V(x) + W , where 
W is a constant that determines the origin of V^\x). Note that because ReW(x) does 
not depend on £ = n/N for n > n*, as shown in Appendix IXj we can put it outside 
the integration with respect to £. We have also used the fact that the integration in 
the second line can be identified with the eigenvalue density p(x), which is shown in 
Appendix El 

Now, we are ready to take the ratio of the partition functions inside and outside the 
cut, and obtain the chemical potential of the instanton. This can be expressed as 

f t . H th f dx(det(x-$) 2 )e- NV ^ . . 

f! . , dx(det(x - $) 2 )e" w ( :r ) ' ^ 

J inside the cut N v ' ' 

Most of the contribution from outside the cut comes from the saddle point of the effective 
potential in the large- limit. In our case, there is a maximum of the effective potential 
at x = 0. This maximum always exists when we consider the case where there are 
two cuts that are symmetric under x — > —x. We consider only the instanton that is 
located at this maximum. Using the method of steepest descent, we obtain the chemical 
potential with an overall factor of O(N ) as 



/ 7T 


r N + r N + 


r N - rjv|| 


N(R'(0) - ±F"(0)) 


47rv / Tv 


r N - r N 





N I dx {2R{x) - V'(x)} 

J a 



V = hl TVTJTTT^ — i ^^exp 

(3.20) 

where a corresponds to Wq and can be chosen as any point in the left one of the two 
cuts. As mentioned above, it indeed determines the origin of V^\x). If there is another 
maximum or minimum of the effective potential, we should take it into account in the 
contribution from outside the cut. Under the double scaling limit, however, there can be 
no contribution from such a saddle point. Hence, it is sufficient to consider the instanton 
at the maximum of x = 0. We will elucidate this point further in the next section. 
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4 Universality of chemical potential 



In this section, we consider the double scaling limit in the case of two cuts. This is 
obtained by taking the limits N —>■ oo and g — > g c with a certain combination of them 
fixed. Here, g is a parameter of the potential V(x) and g c is its critical value. The 
critical point around which the type OB noncritical string is described is one of the 
Gross- Witten phase transitions. It can be found as the critical point of the one-matrix 
model where two cuts, which are the supports of the eigenvalue distribution, become 
closer and are merged. If this critical point is exceeded, we have a one-cut eigenvalue 
distribution. Hence, this critical point distinguishes two phases, namely, the one-cut 
phase and the two-cut phase. 

Before taking up the double scaling limit, we consider the behavior of the matrix 
model near the critical point in the two-cut case. To compute quantities with the method 
of orthogonal polynomials, it is necessary to evaluate the value of r n . This can be done 
by using ()3.4b|) . Because P n (x) is monic, V n>n -i = n. Picking up the coefficient of P n -\ 
in (J3.4b|) . we obtain at the leading order of 1/N 

g£ = F{r,r) = F{r,r). (4.1) 

Here, F(x, y) originates from [<7V'(X)] n In fact, there are two relations derived from 
(13.4b)) . namely, one for even n and one for odd n. If we define F(r,r) for even n, the 
same variable for odd n can be obtained by interchanging r and r, which is nothing 
but F(r,r). At the critical point, the two cuts merge and r and r take the same value 
r c . This means that r = r = r c for n = N at the critical point. Hence, we obtain 
g c = F(r c ,r c ). Expanding (|4.1|) around r c , it can be expressed as 

F(r, r) = g c — A(r c — r ) — A(r c — r) 

- B(r c - r) 2 - B(r c - r) 2 - C(r c - f)(r c - r). (4.2) 

It is necessary for at least one of following conditions to be satisfied so that F(r,r) = 
F(r,r) holds: 

Case 1: r = r 

Case 2: A = A and (r c — r) + (r c — r) = 
Case 3: A = A and B = B. 



12 



The condition that is relevant to the two-cut case is case 2. In this case, if g < g c 
solution of (j4.1|) with real r's exists that can be rewritten as 



r = r c - a^/g c - g£, (4.3a) 
r = r c + ay/g c - g£. (4.3b) 

For g > g c , we cannot take this condition. In this case, the condition r = r will be 
satisfied, which is the case of one cut. Using the definition in 1)4.3)) . a solution of (j4.1|) 
in this case can be expressed as (see Appendix EJ) 

a 2 

r = r e -—(g c -g£). (4.4) 

First, we evaluate the free energy without the contribution from the instanton. In 
the Gross- Witten phase transition comprising the third order, the third derivative of the 
free energy in general has a discontinuity, which is universal in the sense that it is not 
affected by details of the potential. Thus, in order to extract the universal part of the 
free energy, we will compare the free energy of each phase — one cut and two cuts — and 
pick up a discontinuity of (the third derivative of) the free energy, and which is the only 
contribution to the universal part. For the free energy in the two-cut phase, from ()3.3)) 
we obtain 

„ N 2 n 



2 



[ d£(l-£)[logf(£)+logF(£)] (4 . 5a) 
Jo 



N 2 ° f 1 fa 2 ^ 2 



o 

N 2 f a ^ 2 



— J d£(i-£)(-) (9c- g0 (ui.) 



12 \gr. 



(9c - 9) 3 , (4.5c) 



and for the one-cut phase, 

F = 



N 2 /'df (l-£)logr(0 (4-6a) 
Jo 

s -^( i -e(9'(*-«fl 

N 2 fa 2 \ 2 
- -j- J o d£ (1 - I- J (9c - gti) (4.6b) 

--^{gvS^^ (4 ' 6C) 
Here, we have omitted the terms that do not contribute to the discontinuity of the free 
energy. In computation of the free energy in the one-cut phase, we have noticed the 
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following fact: in the one-cut phase, in general there is only one cut constructed from 
N eigenvalues. However, if there were n < n Q = Ng c /g eigenvalues, the cut would split 
into two. Hence, for £ < £ = no/N, we should take the two-cut solution 1)4.8)1 . even in 
the one-cut phase. Comparing (|4.5c)l and (|4.6cjl . we can see that the third derivative of 
the free energy indeed has a discontinuity. 

Second, we consider the contribution from the instanton in a background with a fixed 
number of instantons. In the one-instanton background, an instanton is located at the 
top of the effective potential. The contribution to the free energy from this instanton 
can be obtained by the height of the potential barrier. The effective potential can be 
obtained from the resolvent that is expressed near the critical point as 

1 yl - ~ 1 , 

R(x) = -d x d£ fc (0) (O + fc (0) (O ~C'xVa 2 -x 2 , (4.7a) 

2 Jo L J 

C" = V> a 2 = -(g c -g), (4.7b) 

a z g r c 

where we have dropped the contribution from the potential V(x), which is nonuniversal. 
Then, the height of the potential barrier is 

2N [ dxR{x) = -In— (g c -gf. (4.8) 
J a 3 gr c 

This is the contribution to the free energy from the instanton in the one-instanton 
background. In the n-instanton background, the contribution of the instantons at the 
leading order is multiplied by n. 

Third, we evaluate the chemical potential of the instanton. From (j2.10b|) . (J4.3)) . and 
we obtain 



^ = , r^r ^ ^ exp 



2 » T Oi . .3 
--N—{g c - 9 y 

3 gr t 



c 



(4.9) 



Now, we are ready to consider the double scaling limit. Two limits, iV — > oo and 
g — > g c , are related as 

2 



N = a- 3 [—Y(9c-g) = aH a -> 0. (4.10) 

There is an ambiguity in the definition of t. Here, we define t so that the discontinuity 
of the free energy is described as F"(t) = — \t\/8. Using this definition, the free energy 
in our calculation and that in the string equation coincide. The string equation in this 
definition takes the form of the Painleve II equation, 

th = h 3 -2d?h. (4.11) 
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Here, h(t) is related to the free energy as 

dlF = - l -h 2 + f t = h 2 -Af. (4.12) 

The solution of ()4.11j) in the two-cut phase is given by 

h(t) = Vi+--- (4.13) 

in perturbative expansion. The solution in the one-cut phase is h(t) = 0. Thus, the 
discontinuity of the free energy is consistent with ours. 

In this double scaling limit, the chemical potential of the instanton becomes 



fi = 3 exp 



2 3 

3 



(4.14) 



As shown above, the value does not depend on the details of the potential and hence 
it is a universal quantity. Moreover, it agrees with the nonperturbative effect derived 
from the string equation. This can be seen as follows. If there are two solutions of the 
string equation with the same asymptotic expansion for large t 

h(t) = Vt-~r% + ---, (4.15) 

then the difference between them due to the nonperturbative effect is given by 

hmstit) = Const x r^e-i* 3 . (4.16) 

From this, the nonperturbative effect of the free energy can be obtained as 

3 

F inst = Const x He~i* 5 , (4.17) 

which is in accord with our chemical potential of the instanton. This justifies our 
identification of the instanton in the matrix model as the nonperturbative effect of 
string theory. However, the overall normalization factor cannot be determined from the 
string equation itself. We have determined it from direct computations in the matrix 
model and have shown that it is universal. 



5 Conclusions 

In this paper, we have investigated in full detail the nonperturbative effect in type OB 
string theory, which is defined by taking the double scaling limit of the one-matrix model 
with a two-cut eigenvalue distribution. We have computed the contribution from the 
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instanton to the free energy directly in the matrix model. In the double scaling limit, 
the chemical potential of the instanton does not depend on the details of the potential, 
and is a universal quantity. It takes exactly the same form as the nonperturbative 
effect derived from the string equation. Moreover, by computation via the matrix model 
keeping iV finite, we have fixed the overall factor of the chemical potential, which cannot 
be determined by the string equation itself. In j^, it is shown that in (bosonic) c = 
noncritical string theory, only the closed string, or the string equation, does not describe 
the nonperturbative effect completely, and that the matrix model is more fundamental, 
capturing the nonperturbative effect correctly. Here, we have found that this is also the 
case with type OB, or c = 0, noncritical string theory. 

It is worth noting a crucial difference between c = and c = string theory. In 
the case of c = string theory, the potential of the matrix model is unbounded from 
below. Therefore, the vacuum with one cut is unstable and the imaginary part of the 
free energy obtained in |lj reflects this instability. On the other hand, in the case of 
c = string theory defined by a matrix model with a two-cut eigenvalue distribution, the 
potential is a double-well type and bounded from below. Thus, the matrix integration 
in the definition of the partition function or the free energy as in (j2.1)l is well defined 
and gives a real number. In fact, in the double scaling limit, the integration outside the 
cut in (j2.10b|) becomes 



Therefore, the contribution from the instanton to the free energy is given by an inte- 
gration of a real function, but it diverges. This divergence can be attributed to the 
large- N limit, where some eigenvalues are at the edge of the cut £ = — Moreover, 
in the double scaling limit, some of the eigenvalues are pushed out of the cut ( > —Vt, 
so that the edge of the eigenvalue distribution is smeared. In this sense, the boundary 
at C = — \ft between the outside and inside of the cut becomes subtle, and it becomes 
somewhat artificial in the double scaling limit to divide the integration region into these 
two regions. Because this divergence originates from only a part of the eigenvalues 
that spread into the outside of the cut, it is less divergent compared to the integration 
over the inside of the cut given as the denominator in (|2.10b)l . which is of order N 
as shown in (JHHHJ). Therefore, it may be possible to change the definition of "inside 
the cut" slightly and then renormalize the above divergence into the contribution in 
the denominator in (|3.19j) . In that case, it is important to confirm that the chemical 
potential is still universal irrespective of a slight change in the definition of the inside 
of the cut. However, at least as long as the string coupling constant g s is sufficiently 




(5.1) 
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small, 2 we can take advantage of the saddle point method to evaluate the above integra- 
tion. Thus, the above integration is dominated by the saddle point ( ~ and if we can 
choose the contour of ( as the whole imaginary axis, the integration becomes finite and 
gives an imaginary number as in (j4.14j) . The chemical potential obtained in this way 
is meaningful at least for small g s . In fact, the instanton in the matrix model we have 
considered corresponds to the D-brane known as a ZZ-brane [H]. This can be checked 
by noting that it gives an open boundary to the worldsheet, or more quantitatively, by 
comparing the disk amplitude in the fixed instanton background to that in the ZZ-brane 
background computed in the super Liouville theory as done in jj]. In type OB string 
theory, the ZZ-brane is unstable. The chemical potential that we have computed as 
above is a purely imaginary number. It can be considered to reflect the fact that the 
D-brane under consideration is unstable. Although c = string theory does not have 
the time direction and the instanton does not decay, if we consider the additional time 
direction associated with the energy of the system, the instanton will decay due to its 
instability. The chemical potential of the imaginary number can be considered to show 
this instability. In this sense, the chemical potential is analogous with the energy of a 
statistical system, where there is no time direction, but the energy gives the probability 
that the system is realized in the ensemble. Likewise, because the instanton corresponds 
to the ZZ-brane in c = noncritical string theory, the chemical potential is related to 
the decay rate of the ZZ-brane. Note that it is not necessary for our result to agree 
with the computations in [JJj of the decay rate of the ZZ-brane (D-particle) in c = 1 
noncritical string theory, because the time in the target space is a priori not the same 
as the extra time direction associated with the free energy mentioned above. 

The most important problem remaining unsolved is to identify the boundary condi- 
tions of the string equation. Note that the matrix model from which the string equation 
is derived has no ambiguity. Therefore, it is natural to expect that the matrix model 
specifies some boundary conditions, by which we can determine the nonperturbative 
effect precisely. Because the string equation is a differential equation of the second or- 
der, there are two boundary conditions to be specified. It can be readily seen that the 
asymptotic behavior (j4.15|) needs fine-tuning of one parameter, and therefore there is 
one boundary condition remaining. It is not clear in general whether the contribution 
from the instanton to the free energy, or the nonperturbative effect in the free energy, 
has something to do with this boundary condition. In particular, if g s is sufficiently 
small, the nonperturbative effect is exponentially small and it is impossible in general 
2 g s can be restored on dimensional grounds as i 3 / 2 — > t 3 / 2 /g s . 
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to determine it as an ambiguity to be added to a perturbative series of the free energy, 
which is an asymptotic expansion. However, if we concentrate on the process where each 
term of a perturbative series is completely zero, we can identify the nonperturbative ef- 
fect unambiguously like the tunneling effect. The chemical potential we have computed 
can be regarded as such an example. That is, it is purely imaginary and the free energy 
cannot have an imaginary part perturbatively. We therefore conclude that by defining 
the contour of the integration with respect to the eigenvalue as the imaginary axis, we 
can compute the chemical potential of the instanton, and this will give a boundary con- 
dition for the imaginary part of a complexified string equation. However, we have not 
yet fixed a boundary condition for the real part of the string equation. Note here that 
in order to identify the nonperturbative ambiguity, we should specify a choice of the 
contour in general as mentioned in [T2*] . 

Finally, it would be interesting to apply the computations employed in this paper to 
other matrix models; for example, the two-matrix model. This is left for future studies. 
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Appendix A: Orthogonal polynomials 



In this section, we derive the asymptotic behavior of the orthogonal polynomials both 
inside and outside the cut. This can be deduced using (|3.4|) . From (|3.4ajl . we write the 
orthogonal polynomials in terms of r n . Using the ratio of the orthogonal polynomials 
e kn = P„(x) I P„_i(x) , ()3.4a|) can be expressed as 

x = e kn+1 + r n e~ kn . (A.l) 

Here we have used the fact that s n is identically zero due to the symmetry with respect 
to x — > —x. In the case of two cuts, we should distinguish quantities for even n and odd 
n as follows: 



k n (for even n) J r n (for even n) 

k n (for odd n) [^Vi (for odd ra) 

From (jA.lj) we obtain two relations, one for even n and one for odd n: 



(A.2) 



x = e hn+1 + r n e kn (for even n), (A. 3a) 

x = e kn+1 + r n e~ kn (for odd n). (A.3b) 

In the large- N limit, k n , k n , r n , and r n can be expanded as 

/n+l = /(£) + jfWjj) + ■ ■ ■ > ( A - 4a ) 
/» = /f + + (A.4b) 

Here, /„ represents k n , k n , r n , or r n and /(£) is a continuous function of £ = n/AT 
corresponding to f n in the large- A?" limit. We can see that fi 1 ^ and rfP are identically 
zero from (|3.4bjl . From now on, we will omit the superscript "(0)" for r n and r n . 
Substituting (jA.4j) into ()A.3|) . we obtain 



i = e " + r n e n (for even n), (A. 5a) 

r(o) _ r(o) 

x = e n + r n e n (for odd n), (A.5b) 



at the leading order, and 

= + dfi^ e^ 0> - rJfeWe - *^ (for even n ), (A.6a) 

= f + dfi^ e^ 0> - r n k^e~ k " 0) (for odd n ), (A.6b) 
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at the next-to- leading order. From (|A.5|i . we obtain 



r(o) 



e " 



— | (x 2 - r n + r n ) ± a/ (x 2 - f n - f n ) 2 - 4f n f n } 



(A.7a) 
(A. 7b) 



where the sign of ± is + for x 2 > r n + r n + 2\/r n r n so that e fc ™ ) ~ x as |x| — > oo, and 
— for x 2 < r n + f n — 2y/r n r n not to diverge at x = 0. Moreover, from (|A.6[) . 

+ fcW = --<% log [(x 2 - r„ - F n ) 2 - 4f n F n ] . 
By using these ks, the orthogonal polynomials can thus be expressed as 



(A.8) 



P n (x) = exp 



22 kn 



m=l 



(A.9) 



Because we need to perform computations up to the next-to-leading order of l/N, we 
use the Euler-Mclaurin summation formula to obtain 



Pjx) 



exp 



2 Jo+1, 2 



P„(x) = exp 



N 



N 



W T AT 



' N ~ 



(for even n) 
(for odd n). 



(A. 10a) 
(A.lOb) 



Substituting (TOT) and (TOl) into (lA~Tol) . we find 



P n {x) 

W (x,C) 
M 



x 



Here, kn' is for even n and /cn Uj for odd n. 

If x is in the oscillatory region of P n (x), where (x 2 — r n — r n ) 2 — 4/r n r n < 0, k n is 
no longer real, but P n (x) should be real. We note here that P n (x) is a solution of a set 
of linear differential equations ()3.4|) and we have obtained its solution in a classically 
allowed region, (jA.llb|) via the WKB approximation. Therefore, in order to find P n in 
the oscillatory region which is an analog of a forbidden region, we invoke the continuity 
condition jl]. That is, it is sufficient to take a linear combination of two complex 
solutions of ()3.4a|) to obtain the real solution, which is nothing but the real part of the 



V A ' n 1 n) ^' n' n 

d£' (fc (0) (f) +k (0) (C 



(0) 



exp 



N 

T 



(Alia) 
(A.llb) 
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naive expression in (jA.llb|) with an appropriate phase factor. We thus obtain P n (x) in 
the oscillatory region as 



Pn{x) 



X TV 



^" n' n 



(x 2 - r„ - rj 2 



x exp 

Here, we have used 



N / n 

-rReWnlx,- 



sm 



Im(^W (x^) + h£A+e}. (A.12) 



Re^°) = ilogr n , (A.13) 



in the oscillatory region. The additional factor 2 comes from a relative normalization to 
the solution for the allowed region (x 2 —r n — r n ) 2 —Ar n r n > in the continuity condition, 
and 9 is a constant phase factor. The quartic root is defined as real and positive, and 
the overall sign is included in 9. 

If we consider the normalized orthogonal polynomials ip n (x) = P n (x)/\^h^, This can 
be expressed as 



X 2 Vn 



kr n r n - (x 2 -r n — r n ) 2 
x exp 



WAx 



n 



W [x 



f ReW r (*, £)] sin [im (^W (x, £) + + 9 

d£' (log f(£')+ log f (£'))• 



n 
N 



, (A. 14a) 
(A. 14b) 



For n > n*(x), where n*(x) is defined by (x 2 — r Ut — r n J 2 — 4r n<t r n „ = 0, n satisfies 
n > if x is in the oscillatory region of P n (x). From ()A.13|) . for n > n*, 
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ReH/ r (x,£ = - 



d£ ( £ (0) (0 + A: (0) (0 - 5 logr n - i logr n 



(A.15) 



hence ReH^x, £) does not depend on £. Thus, we can set W r (a;, £) = W r (x, 1) for 
£ > n*/N. As seen in (jA.12|) . the frequency of the oscillation of P n (x) is of order N. 
Therefore, if we take the average of this oscillation for P 2 (x), it simply gives 1/2. We 
thus obtain 



y/A? n r n - (x 2 -r n - r n )'- 



:exp[NReK r (x,l)]. 



(A.16) 



Here, the branch of the square root is chosen to be real and positive. This expression is 
convenient for calculating because of the independence of £ at order N. 
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Appendix B: Effective potential and resolvent 



In this section, we clarify the relation between the effective potential and the resolvent. 
We show that the first derivative of the effective potential can be identified with the 
resolvent. The coefficient of the resolvent in the double scaling limit can be determined 
by this relation. 

The relation between the effective potential and resolvent can be derived as follows. 
Because Pn(x) = (det(x — $)), we have 



„2 X 4 



fB.r 



for x outside the cut. Differentiating both sides with respect to x, we obtain 

jf tr ^3$ det ( x -^) = < det (* - $ )) (\ d ' K ( x ) + °(^)) ■ ( B -2) 

In the large-iV limit, the expectation value on the left-hand side factorizes. Hence, we 
have 

= \d x K(x) = \d,f it (fc<°>(£) + k'°>(0) . (B.3b) 

The resolvent inside the cut can be obtained by analytic continuation, and it is consistent 
with the choice of the double sign in the definition of fc(°)(f) and fc(°)(f) given m (I3l6bl 
and (jSHEcJ). We can evaluate the integration in ()3.f9|) by using this relation. In the 

case of two cuts, we have 

d x (fc<°>(0 + * (0) (0) = 7( 2 5 = , 2 =^=, (B.4) 

v J \f{x-r n -r n y-Ar n r n 

which is again understood as an analytic function. Thus, the eigenvalue density can be 
obtained as 

p( x ) = — /l-tT _)--?-/ tr -\ (B.5a) 

1 Re f d£—j= X (B.5b) 
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Here, the square root is again defined as an analytic function. The sign of the square 
root is positive for positive x. Note that p(x) is positive even for negative x because 
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the square root also becomes negative there. Thus, the integration in (|3.19|) gives the 
eigenvalue density. 

The relation between the effective potential and the resolvent ()B.3j) can be checked 
by a direct computation near the critical point, where it can be expressed as 

d x (k {0) (0 + k {0) (0) ^ ; X (B.6a) 

a 2 

a(0 2 ^—(9c- 90 b 2 (0 ~ 4r c . (B.6b) 

r c 

Here, a = a(l) and b = b(l) are identified with the endpoints of the cut. Integrating 
this equation with respect to x, the effective potential near the critical point can be 
obtained as 

\v£°\x) = -\d x J d^ [k(0 + k(0] * -C'xV^^, (B.7a) 



(" ^ a 2 = -(g c -g), (B.7b) 

a 2 g r c 



and can be identified with the resolvent. 



Appendix C: Relation between r in one-cut phase 

and in two-cut phase 

There is a relation between the behavior of r in the one-cut phase and in the two-cut 
phase. Here, we use the notations in (|4.2j) 

g £ = F(r, r) = g c — A(r c - r ) - A(r c - r) 

- B(r c - r) 2 - B(r c -r) 2 - C(r c - f)(r c - r), JOJ) 

and from (|4.1jl we consider the case where A = A = A. In the two-cut phase, r and r 
satisfy the condition (r c — r) + (r c — r) =0. From these conditions, r near the critical 
point in the two-cut phase is expressed as 

r = r c ± [B + B-CY~ 2 ^g c -gi. (C.l) 

The double sign indicates + for r and — for r, because r vanishes for C, = 0. In the 
one-cut phase, r = r = r , and then r near the critical point is 

r = r c -±-Ag c -gt). (C.2) 
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A, B, B, and C can be expressed in terms of F(r,r) as 



A 
B 
B 

C = - 



dF(r, r) 



dF{j, r) 



dr 

1 d 2 F(r,7) 



dr 



2 dr 2 
1 d 2 F(r,7) 



2 dT 2 
d 2 F(r, r) 



drdr 



(C.3a) 
(C.3b) 
(C.3c) 
(G.3d) 



Here, we use the fact that F(r,r) originates from [gV(X)] nn _ 1 , and that it can be 
expressed in terms of V(x). Because the potential we consider is even, V'(x) has only 
the odd power of x and can be expressed as gV'(x) = xll(x 2 ). Using this, F(r,r) is 
expressed as 



F(r, r) 



dz[z + -)u((z + -)(z + -)), 



(C.4) 



for even n, and whereas odd n we should use Fir, r). The condition A = A is not trivial, 
but comes from (|4.1jl . Using ()C.4j) . (|4.1j) is expressed as 



f-r) (b —U{[ z+'- \ [z + - 





To satisfy this condition with f ^ r, we need 

dz 



(C.5) 



= G(r,r) 



U((z + 



z + 



(C.6) 



From (£22), (KI4T1 . and (KlBTl . we can see 



B+B-C = — 
2r r 



(C.7) 



This is the relation between r in the one-cut phase and in the two-cut phase. 
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